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Abstract 

We show that the behaviour of the outgoing radial null geodesic 
congruence on the boundary of the trapped region (suitably defined 
as a four dimensional region) is related to the property of nakedness 
in spherical dust collapse. The argument involves a conformal trans- 
formation which justifies the difference in the Penrose diagrams in the 
naked and covered dust collapse scenarios. 
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1 Introduction 



Consider a cloud of matter (regular initial Cauchy data) collapsing in- 
definitely under its own gravity. A singularity eventually develops in 
the spacetime and it is indicated by the divergence of the Kretschmann 
scalar. In the advanced stages of collapse trapped regions are formed 
[1], [2] and there exists a null ray which marginally escapes to in- 
finity (event horizon). It is not clear whether the singular boundary 
is entirely surrounded by the trapped region. In other words, it is 
not known if a portion of the singular boundary is exposed in the 
untrapped region and non-spacelike geodesies can emanate from it 
(naked singularities). In fact, exact solutions to Einstein's field equa- 
tions with certain kind of source terms are known to exhibit both 
naked and covered singularities depending upon the sort of regular 
initial data chosen (See [3] for details). 

The complexity of the general problem lies in the fact that the 
Cauchy initial value problem for the Einstein field equations with 
sources is less tractable. The systematics available about the initial 
value problem is far too small for any implication for questions like 
formation of naked singularities. For instance, even the well-posedness 
of the problem is not self evident and has to be proved independently 
for different types of sources [4], [5], [6]. It would be indeed difficult, 
for example, to find or even expect a conserved or monotonically be- 
having function of the Cauchy surface with respect to its evolution, 
which could be expected to provide insight into the process of creation 
of naked singularities. 

As a result of this difficulty, a large number of investigations that 
have been carried out, have been concerning certain exact solutions 
or numerical simulations ([3] and references cited in them). Issues 
like strength of the singularities, genericity, behaviour with respect to 
change of source etc. have been studied in examples like dust, null 
dust, perfect and imperfect fluids and scalar fields. However, they do 
not suggest any typical geometrical feature which could be expected 
to arise before a naked singularity forms ( In case of a singularity the 
singularity theorems make use of a typical geometrical feature viz. 
trapped surfaces to prove its existence). The lack of indication of 
existence of such a feature is evident in the fact that one is forced 
to check for the existence of naked singularities in a direct manner 
whenever required. To be precise, one simply checks if non space- 
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like geodesies emerge from the singular boundary using differential 
geometry and the form of the metric in the example. 

This paper is a first step towards an indirect criterion. Preferably, 
the criterion should be applicable away from the singular boundary. 
That is a non-local problem and given the difficulties of Cauchy evo- 
lution, there is no indication available for what the criterion could be. 
It would be perhaps appropriate in such a situation to examine regions 
near the singularity. There are also parts of the singularity from which 
geodesies cannot escape. Some criterion applicable to such a portion 
would also be significant. It would indicate that the information about 
the exposure of a part of the boundary is contained elsewhere on the 
boundary. 

This work is the first part in which the example of spherical dust 
collapse is considered. This example is well studied and the initial 
data which leads to naked singularities is known. 

2 Spherical Dust Collapse 

We illustrate the main features of spherical dust collapse in radial co- 
ordinates (figures 1 and 2) and in Penrose diagrams (figures 3 and 
4). 

Point O is defined as the point where the apparent horizon 2 meets 
the singularity. In figure 3, there is no portion of the boundary be- 
yond O that is exposed. In figure 4, however, there is a null portion. 
In this paper, we examine if O would yield the information about the 
existence of an exposed portion. O is a covered point of the singular 
boundary (with or without a naked portion). Mathematically, one 
works with points on the apparent horizon in the approach to the 
concerned point O. This could be looked upon as a property of the 
marginal trapped regions which constitute the apparent horizon. One 
is therefore working with outgoing null congruences with zero expan- 

2 In this work we define a trapped region as the union of all possible closed trapped 
surfaces (See [5] or [4] for the definition of a closed trapped surface). This is in general a 
four dimensional region. The boundary of this region within spacetime is referred to as 
'apparent horizon' in this work. However, according to the usual definition (See [5] and 
[4]) , a section of this boundary by a Cauchy surface would be called an apparent horizon. 
We make this change for lack of established terminology for the ideas we use (which are 
independent of any Cauchy surface). See also figure 60 of [5] and [10]. 
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Figure 1: Collapse of spherical dust leading to a covered singularity 



Singularity 




Figure 2: Collapse of spherical dust leading to a naked singularity 

sion. If a singularity meets such a region, the point in the strict sense 
will not have any non-spacelike geodesies emerging from itself and will 
therefore be covered. 

We proceed as follows. Consider a congruence of outgoing null 



4 



Singularity 
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Figure 3: Penrose Carter diagram for collapse of spherical dust leading to a 
covered singularity 

geodesies. Let us parametrize each geodesic using an affine parameter. 
We investigate the tangent vectors £ to geodesies in the congruence 
at the apparent horizon (£ M = dx^/ds where s is an affine parameter 
along the null geodesic of the congruence) . 

O is a part of the singularity which is marginally trapped (expan- 
sion of any outgoing congruence will be zero here) and so there will be 
no outgoing geodesic emerging from O. Any metric which confirms to 
that is referred to in this paper as having the 'correct' causal structure 
at the boundary. 

In the limit of approach to O, the vector £ is expected to vanish 
if the metric used depicts the appropriate causal structure at O. O 
is a part of the singularity which is marginally trapped (expansion 
of any outgoing congruence will be zero here) there will be no out- 
going geodesic emerging from O. Any metric which confirms to that 
is referred to in this paper as having the 'correct' causal structure at 
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Figure 4: Penrose- Carter diagram for collapse of spherical dust leading to a 
(locally) naked singularity 

the boundary. We examine this in the case of spherical dust collapse 
where the metric is evolved with regular initital Cauchy data. It is 
found that the condition is met at O for the collapse leading to covered 
singularities. However, the tangent vector stays non-zero in the limit 
for the naked case. This is remedied by performing a conformal trans- 
formation which diverges at O. The appropriate causal structure is 
then restored. Thus, in case of spherical dust collapse, this behaviour 
of the tangent vector is related to the property of nakedness of the sin- 
gularity. Could this be true in general? It would mean that a diverging 
conformal tansformation is necessitated at the singular boundary in 
the naked case as against the covered case. Some theorems in the case 
of general collapse are provable and work is in progress. 

Thus, it is suggested that in a general collapse, the tangent vector 
field of the congruence on the apparent horizon (precisely, its limiting 
behaviour at the trapped part of the singularity) contains the informa- 
tion about whether the singular boundary has a naked portion. 

The plan of the rest of the paper is as follows. We wish to in- 
vestigate the tangent vector on the apparent horizon. We choose the 
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radial null geodesic congruence in spherical symmetry. The general ex- 
pression for tangent vector is not available in a closed analytic form. 
However, one can calculate it in a simpler special case of dust, the self 
similar model, and subsequently show how to generalize the results. 
The next section describes the self similar dust model, where we dis- 
cuss the tangent vector and demonstrate the connection with naked- 
ness. In the next section we show that the results can be extended 
to the general dust case. We then turn to the mentioned conformal 
transformation in the next section. It leads to the Penrose diagram for 
the naked case from the usual metric of Tolman and Bondi in spher- 
ical co-ordinates. We show that the conformal transformation should 
diverge at the singular boundary. This results from the fact that the 
tangent vector ought to vanish at O in the correct causal metric. 



3 Self Similar Dust Model 

The collapse of a spherical cloud of pressureless fluid is given by the 
following metric [7]. 

ds 2 = dt 2 ^vrrdr 2 - R 2 cM 2 (1) 

l + /(r) 

where 

t and r are the co-moving time and radial co-ordinates respectively. 
R(t, r) is called the 'area radius' and a closed expression for this quan- 
tity which results in the dust case has enabled substantial progress in 
understanding the model. 

Two free functions arise viz. F{r) , called the mass function since 
it is the total mass to the interior of a shell of radius r and the total 
energy function f(r) which is called so because of the constraint below 
which resembles a relation between kinetic and gravitational potential 
energies of a shell. 

R2 = ^ + f ( 2 ) 

The source energy momentum tensor is diag[p(t, r), 0, 0, 0]. 
The solution for R mentioned above is 

t- H r) = -^G{=f) (3) 
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where a singularity boundary is formed at t = to(r). The function 
G is defined as follows. 



G(y) 



™^_v^ l> y >0, 

§, _ y = o, 

-axcsinVEy _ VIEv >y>-OO. 



(_„)3/2 - y 

The central shell focussing singularity, which is the limit as r — > 
along this locus is of interest and turns out to be naked for some initial 
data. 

The self-similar model is the one in which F(r) = Xr where A is 
a constant (which decides if the central singularity will be naked or 
not) and f(r) = 0. 3 

We choose the scaling to(r) = r. A self similar co-ordinate z = t/r 
is introduced. We note the expressions for R and R' which will be 
useful in the subsequent analysis 

R = r\- 2 / 3 (3/2{z-l)) 2/3 (4) 

*-(25f(^) 

We cast the metric into double null co-ordinates. It is not difficult 
to show that 

ds 2 = r 2 (z 2 - i?' 2 ) du dv (6) 

where 

dr dz ._. 

dr dz 
dv = — + — — (8 

r z + R'{z) v ' 

The double null form (ds 2 = C 2 (u,v)dudv) turns out to be useful 
when affine parameters along null geodesies are to be calculated. For 
instance, along an outgoing radial null geodesic (du = 0), the affine 
parameter is J u=consta nt C 2 dv up to a multiplicative and an additive 
constant. 

Now let us turn to calculating the tangent vector to the outgoing 
null radial geodesic congruence, which is our primary interest. 



3 The model is referred to as self similar since there exists a nomothetic Killing field 



' dt ^ 1 dr 
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Assume the vector to be of the form 



Z= (Q(t,r),Q(t,r)y/T+f /X,0,0) (9) 

where Q is obtained from the geodesic equation which £ has to 
satisfy. That constraint turns out to be 



QQ + QQ'VT+f/R' + Q 2 R'/2R' = (10) 

We have provided the expressions for the most general dust case 
here. One may read off the expressions for the self similar case by 
setting / to zero and using equation (5) for R'. 

The equation above takes the form 

1/Q = I ^-dk + A(u) (11) 

Ju=constant 

where k is an affine parameter along outgoing radial null geodesies 
and u is the retarded null co-ordinate. A is an arbitrary function of u 
resulting because of the partial integration. 
As indicated earlier, 

dk = r 2 [z 2 - R' 2 ) dv (12) 



keeping u fixed. 

Using this in equation (11) we obtain 



VQ= I 

J u- 



r{z- 3f 
3 z 2 -l 



2A/3 



I u=constant 

or using the fact that du = from equation (7) 



1/3 

dz (13) 



u=constant 3 Z 2 1 



2A/3 



z-1 



1/3 



b-R'(b) 



dz (14) 



The integral over z is to be evaluated from r = to the apparent 
horizon, where we shall be interested in evaluating the tangent vector. 



4 The expression for r thus obtained in terms of z has a positive multiplicative constant. 
Whether r vanishes or not( and later whether 1/Q vanishes or is non-zero) is of importance 
in the calculations. So setting this constant to 1, which is done here, will not affect the 
result. 
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The latter can be shown to be the curve R = F and turns out to be 
the locus z = 1 — 2 A/3. 

The integral being over z, it is important to know which value of 
z along the outgoing null curve yields r = 0, the lower limit of the 
integral. This issue as it is shown further leads to the difference in the 
behaviour of Q in the naked and covered cases. 

Consider then the equation 

C db 



which we examine for r = 0. That will happen when the integral in 
the square bracket diverges positively. Two cases can be immediately 
seen to arise. 

Case(i) b — R'(b) =0 has no real root. 
The integrand therefore does not diverge anywhere and also remains 
positive (or entirely negative) all over the real line. It can be checked 
that b— R'(b) > for any one real b which would be sufficient to claim 
that the integrand is positive. Also, b—R'(b) is bounded since b is to be 
limited to the nonsingular region z < 1 (z = 1 is the singularity curve 
itself). So, in order that the integral diverge, the range of integration 
should be infinite. We have chosen to limit the final point to the 
apparent horizon z = 1 — 2 A/3 and hence the initial point must be 
z = — oo. 

In fact a shorter intuitive argument is possible. It is known from 
the Tolman Bondi dust model that central singularity forms at (to(0),0) 
If one assumes the Penrose diagram for the covered case (which is 
indeed what this case turns out to be), the null rays crossing the ap- 
parent horizon begin at the centre at t < to which makes z = — oo 
there. 

Case(ii) b — R'(b) = has at least one real root. 

In this case, the equation implies that r = at the value of z for 
which the integral in the exponent diverges. The range of integration 
for equation (14) would then be limited at the lower end by that value 
of z. This will be the root (in fact the one closest to 1 — 2A/3, the 
apparent horizon). 

That this is so is seen as follows. Equation (15) can then be re- 
cast using an expansion for the integrand in the exponent as outlined 
below. 
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Expanding R'(b) using the Taylor series about the root (called z_) 
it can be shown that the leading order behaviour of r is as follows 



r = (z 



z_)i-« + 0(z - z-) 



(16) 



where 



dR'(b) 
db 



b=z- 



(17) 



(It can be easily checked that a < 0) 
At z = Z-, therefore, r vanishes. 

Thus in conclusion of this analysis, we note that the lower limit of 
integral for l/Q differs. It is — oo when R'(b) — b = can never have 
a real solution and is the root (closest to apparent horizon) when a 
solution exists. 

This observation plays the key role in further analysis. Making 
note of this consider equation (14). Analyzing the various factors in 
the integrand one finds that the integrand would diverge if z = — 1 ( 
z 7^ 1 since we are not on the singular boundary). 

z will take the value —1 in case i. In case ii, the following takes 
place. Consider b — R' using equation (5). It is easy to see that 
b — R' < for all b < 0. So, the root 6_ cannot be negative. Hence it 
is certainly greater than —1. Thus z cannot take the value —1 in case 
ii in the integral for l/Q. 

Thus the integrand diverges as l/(z + 1) in case i and is finite in 
case ii. 

Expanding the rest of the integrand factor in a Taylor series about 
z = — 1, one can easily check that the integral diverges logarithmically 
in case i while staying finite in case ii. 

Thus, Q vanishes in case i and stays non zero (and finite) in case 

ii. 

Returning now to equation (9) , we can now see that £ behaves 
in different ways in case i and case ii on the apparent horizon, in 
particular as one approaches the point O on the Penrose diagrams 
shown (figures 3 and 4) 5 . It can be checked that the factor ^1 + f/R' 
tends to is a non-zero finite quantity on the apparent horizon. 

So, £ vanishes in case i and tends to a non-zero quantity in case ii. 

5 Figure 4 is a diagram of a locally naked singularity. The self similar cloud which we 
examine here turns out to be globally naked. However, the structure near O is the same 
as any locally naked case and figure 4 can be used. 
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Thinking of £ as 'flux density' of the congruence, we may interprete 
that the congruence tends to cluster in case ii as against case i. 

From previous analysis of naked singularities (self similar cases) 
using analysis for emergence of geodesies (roots analysis), it can be 
checked that case i corresponds to the covered case and case ii corre- 
sponds to the naked singular metric. 

4 Extension to the general dust case 

In the general dust case, the equation (10) yields no closed analytic 
solution which would have clearly been useful. However, we note that 
we are interested only in the behaviour of Q in the limit of approach 
to point O on the apparent horizon. 

To this end the following observation plays an important role. It 
is shown that given a dust solution, one can construct a modified dust 
solution (modified distribution) which in a suitable limit approaches 
the given dust solution [8] . The key result that makes this construction 
useful is that it is proved that naked modified distributions reproduce 
naked dust solutions given and covered modified distributions repro- 
duce covered ones. One can then work with the modified distribution 
for the given dust solution and take the limit which preserves naked 
or covered nature. We outline the construction in [8] below 

a) Marginally bound case (/ = 0) 

Imagine a shell of radius r c in the given Tolman Bondi dust model. 
Replace the interior of the shell by a self similar dust metric, matching 
the first and second fundamental forms at the interface r = r c . It 
can be shown that this restricts the self similarity parameter A which 
appears in the mass function. This specifies the self similar solution 
completely. Now taking the limit as r c tends to zero, one can show [8] 
that the matching constraint does imply that the interior self similar 
solution stays naked in the limit if the original dust solution was naked 
and likewise in the covered case. 

b) Non Marginally bound case (/ / 0) 

The construction is similar in this case except for an additional 
interface. Two shells, r c l and r c 2 ( say r c l < r c 2 ) are now considered. 
To the interior of r c l, we replace by a self similar metric. Between 
the two shells, we replace with a dust portion having / so behaved 
that it increases smoothly from zero at r c \ to f(r c 2) of the original 
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dust metric. The F function for this extra portion of dust however is 
the same as that of the original dust metric. We match the first and 
second fundamental forms at each of the interfaces. As before, this 
can be shown to constrain the interior self similar solution uniquely 
given r c l and the original dust solution. Again, the property of being 
naked or covered is preserved in the limit ( r c 2 — > 0) like the previous 
case [8]. 

We now consider Q in the modified distribution for any given dust 
solution. In the self similar part of the latter, results of the previous 
section apply. Since the congruence of outgoing geodesies is smooth, so 
is Q. This makes Q continuous across the interface/s in the modified 
distribution. Now imagine the given dust solution as the limiting case 
of the modified distribution. In the limit of approach to point O on 
the apparent horizon, one has to evaluate Q in the self similar part. 
Because of continuity of Q, the same behaviour will continue to hold 
in the limit of the interface/s tending to zero when the original dust 
solution is reproduced. Making use of the fact that the property of 
being naked or covered is preserved in this limit, one concludes that 
the behaviour of Q (and hence that of £)in the self similar naked and 
covered cases continues to hold in the general dust scenario as well. 

5 Conformal transformation and Pen- 
rose diagram 

The issue about the metric being causally appropriate is related to 
the behaviour of £ by the simple argument below. 

The point O being the intersection of the apparent horizon and the 
singular boundary implies that any outgoing null congruence ought to 
have zero expansion in the limit of approach to O. So O ought not 
to have any outgoing geodesies emanating from itself. Hence £ ought 
to vanish if this causal property is correctly reflected in the metric. 
If the metric is not the correct one, then one performs a conformal 
transformation which diverges at O. This makes the £ vanish and 
yields the causally correct metric. 

We now show that the under a conformal transformation which 
diverges in the limit of O, £ which tends to a non-zero limit transforms 
to a vector field which vanishes in the limit. 

Recall that we defined £ for any geodesic congruence using an affine 
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parametrization. Under conformal transformations, affine parameters 
along null geodesies change ( unlike timelike geodesies which do not 
remain geodesic curves, null geodesies do stay so provided the affine 
parameter changes appropriately). Infinitesimal parameter ds trans- 
forms to ^(x^ds [4], where fi 2 is a conformal transformation. Thus 
it is obvious that £ M = dx^/ds if finite and non- vanishing in the limit 
will vanish under Q 2 transformation provided the latter diverges there. 

Thus we find that in the dust case, one requires a conformal trans- 
formation which diverges on the apparent horizon in the limit of ap- 
proach to the singularity in the naked case as against the covered case 
where the usual metric given in equation (1) is appropriate to describe 
the singularity structure. 6 

This justifies the difference in the structure of the singular bound- 
ary near O in figures 3 and 4. In the naked dust case, when one uses 
the metric used in the previous sections (figure 2), it can be shown 
that several radial null geodesies appear to emerge from the central 
singularity with the same tangent vector [9]. Thus the limiting be- 
haviour is not correct for the set of outgoing radial null geodesies to 
form a congruence, i.e. the condition that there exists only one curve 
of the set passing through any point is violated in the limit of ap- 
proach to the central singularity. The causally approriate structure of 
the central singularity in the naked case is therefore not a point, but a 
one dimensional curve. The diverging conformal transformation does 
precisely that leading from figure 2 to figure 4. The central singular- 
ity being transformed to the one dimensional naked singular portion 
shown) Thus this naked point on the central singularity is depicted 
as a curve, to obtain the Penrose diagram. No such divergence is re- 
quired, in the covered case while going from figure 1 to figure 3 as is 
also implied by the tangent vector vanishing at O in the usual Tolman 
Bondi metric of figure 1. 

6 It is evident from the definitions that the difference in the behaviour of tangent vector 
is due to difference in the metrics, one being conformally related to the other. This raises 
the question of the tangent vector being an appropriate characterization of the causal 
structure. It certainly is not appropriate at any event within the spacetime but in the 
limit of approach of the boundary its behaviour indicates if a metric with the correct 
limiting causal structure has been employed in its calculation. 
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6 Arbitrariness in the definition of £ 



As defined earlier, £ is a tangent vector, defined using an affine param- 
eter, to a congruence of null geodesies the congruence being outgoing 
to a marginal trapped surface in the apparent horizon. Thus, the £ 
field is not unique. The arbitrariness is due to two reasons. Firstly, 
the congruence is not uniquely defined. Secondly, the affine parameter 
used in the definition is arbitrary upto a multiplicative and additive 
constant. 

In spherical symmetry, the first difficulty is remedied by choosing 
a null congruence which is radial on the apparent horizon. This choice 
is not possible in general. It will be shown in the sequel to this paper 
that vanishing (or nonvanishing) of the vector field at O is independent 
of the choice of the null congruence. 

The second difficulty brings in an arbitrariness of a scalar multiple 
of the vector field. This scalar is constant along any null geodesic of 
the congruence. In the approach to the boundary point O, we claim 
that this will not affect the vanishing (or non vanishing) of the field. 
The argument is as follows. 

Consider a vector field £1 which does not vanish at any boundary 
point P. Let (j>(x) be the value of the scalar at any point x on the man- 
ifold which multiplies this vector field (because of the multiplicative 
arbirtariness of the affine parameter). <p is constant along any null 
geodesic of the congruence. We assume that P has a neighbourhood 
which intersects the physical spacetime (spacetime without boundary) 
in a convex normal neighbourhood. So, £1 not vanishing at P would 
imply the existence of a null geodesic beginning from P 7 . Now, if <fi 
vanishes at P then it should vanish all along the portion of the null 
geodesic in the convex normal neighbourhood. This would make the 
multiplicative scaling of the affine parameter zero all along the null 
geodesic, which is impossible. So <\> cannot vanish in the approach to 
P. This completes the argument. 

7 Summary and Conclusion 

The tangent vector field to an affinely parametrized null geodesic con- 
gruence is examined for behaviour on the apparent horizon (defined 

7 This is similar to Thm 8.2.1 of [4] 
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as the boundary of the four dimensional trapped region as in the first 
section) in the approach to the singularity (point O) in the dust col- 
lapse model. There is a correlation with the property of nakedness 
with this behaviour. Demanding that the vector vanishes at the cov- 
ered point O forces the divergence of the conformal transformation at 
O which leads to the Penrose diagram for the naked scenario. Since 
the vector vanishes in the covered case, there is no such divergence 
and hence the Penrose diagrams in the two cases differ. 

Thus, we have shown that the information about whether the sin- 
gularity formed in collapse is naked is contained at the intersection of 
the boundary of the (four dimensional) trapped region and singular 
boundary in the spherical dust case. 

The procedure of checking if an appropriate conformal transfor- 
mation is necessary does not involve checking for emergence of causal 
curves from the singularity. So far, the latter has been the only method 
for checking if a singularity formed in collapse is naked. The work pre- 
sented in this paper suggests an alternative method and demonstrates 
its validity in spherical dust collapse. 
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